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Abstract 

We investigate the one- to two-dimensional zigzag transition in clusters consisting of a small num- 
ber of particles interacting through a Yukawa (Debye) potential and confined in a two-dimensional 
biharmonic potential well. Dusty (complex) plasma clusters with n < 19 monodisperse particles 
are characterized experimentally for two different confining wells. The well anisotropy is accurately 
measured, and the Debye shielding parameter is determined from the longitudinal breathing fre- 
quency. Debye shielding is shown to be important. A model for this system is used to predict 
equilibrium particle configurations. The experiment and model exhibit excellent agreement. The 
critical value of n for the zigzag transition is found to be less than that predicted for an unshielded 
Coulomb interaction. The zigzag transition is shown to behave as a continuous phase transition 
from a one-dimensional to a two-dimensional state, where the state variables are the number of 
particles, the well anisotropy and the Debye shielding parameter. A universal critical exponent for 
the zigzag transition is identified for transitions caused by varying the Debye shielding parameter. 
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I. INTRODUCTION 



Consider a strongly-coupled, two-dimensional (2D) system of n particles with identical 
mass m and charge q . A confining potential well is required to balance the repulsive 
interparticle force and create a stable configuration. For almost any 2D potential well 
expanded around its minimum, the lowest order terms in a particle's potential energy are 
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where Uq is a constant. Consequently, we can approximate the confining potential energy as 

111 1 

U (x, y) = -k x x 2 + -k y y 2 = -muj 2 x x 2 + -mu) 2 y 2 , (2) 



where k x and k y are force constants, and Uq x and ujq v are single-particle (center-of-mass) 
oscillation frequencies in the x and y directions, respectively. In 2D, the biharmonic well 
[Eq. ([2])] gives the general confining potential energy when higher order terms are negligible. 

When charged particles are in free space, they interact through an unshielded Coulomb 
potential l|, |2|, y| . However, if the particles are in a dielectric, then the Coulomb interaction 
is shielded by the medium's dielectric response, and particles interact through a Yukawa 
potential (i.e., a shielded Coulomb or Debye potential), 



V(r) 



47re r 



(3) 



where r is the separation distance, and A is the Debye length. We call a system of particles 
confined to two dimensions and interacting through a Yukawa potential a "2D Yukawa 
system". If n is small, then the system is a 'Yukawa cluster". In 2D Yukawa systems the 
finite Debye length allows the particle-particle interaction length to be varied from long 
range to short range 4J], affecting both the system's static and dynamic properties. 



The case of 2D Yukawa systems in isotropic (uq x 



plored extensively, both theoretically 
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potential wells 



7] and experimentally 



ias been ex- 
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12| . For 



isotropic wells, large-n systems form a circular disk where the interior of the disk has a 



triangular lattice 



2|. For clusters (small-n systems), different shell configurations become 



stable as the interaction is tuned from long range to short range. When the potential well 
is anisotropic, qualitatively new types of configurations can occur 0,B,Q. if tt 
weakly anisotropic then clusters are elliptical and have well-defined shell structures 
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On the other hand, when the well is highl y a nisotropic, the particle configuration is a one- 
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dimensional (ID) straight line 
a zigzag transition. Zigzag co nfig urations may become elliptical and then circular through 



19] . A ID cluster becomes a 2D cluster through 



further structural transitions 



Dusty (complex) plasma should be an ideal experimental system for studying the zigzag 
transition in 2D Yukawa clusters. In laboratory dusty plasmas, monodisperse dust particles 
interacting through a Yukawa potential [20| are confined near the sheath edge above a 
horizontal electrode to form a 2D system. A rectangular confining structure placed on top 



of the electrode can produce a biharmonic potential well 



15 



16 



Melzer 



2l|. 



21| experimentally observed zigzag transitions in dusty plasmas confined in a 



radio frequency (rf) discharge as a function of particle number n and neutral gas pressure, 
and attempted to infer the well anisotropy and Debye shielding parameter using a static 
analysis of the cluster configuration together with a comparison to unshielded Coulomb 
theory. He concluded that the measured cluster properties, including the critical value of n 
for the zigzag transition, were not inconsistent with the physics of an unshielded Coulomb 
interaction (i.e., A — > oo). 

In this paper, we study Yukawa clusters in one- and two-dimensional configurations and 
the transition between these configurations. Dusty plasma experiments are performed as a 
function of particle number for two rectangular confining wells, giving two values of the well 
anisotropy. This work extends previous experiments 15| on 2D Yukawa clusters in weakly 
anisotropic wells. We directly measure the well anisotropy and Debye shielding parameter 
[ijj]. We find that Debye shielding is important, i.e., our results are not consistent with 
physics in the unshielded Coulomb regime. Using the measured cluster parameters, we 
compute predicted equilibrium configurations from the model of Sec. [TTJ The predicted 
and measured configurations exhibit excellent agreement. The critical value of n for the 
zigzag transition is found to be the same in both experiment and theory, and to be less than 
that predicted for an unshielded Coulomb interaction. Even though n is small, the zigzag 
transition is shown to behave as a 1D-2D continuous phase transition and a universal critical 
exponent is identified. 
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II. MODEL 



Two-dimensional Yukawa clusters can be modeled as a strongly-coupled system of n iden- 
tical particles with charge q and mass m at positions {xi,yi} interacting through a Yukawa 
potential [Eq. Q] with Debye length A. The particles are confined in a 2D biharmonic 
well [Eq. ()2])] where u>o x and uo y are oscillation frequencies for the x (longitudinal) and y 



(transverse) directions, respectively. The separation distance between partic 
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es i and j is 



15| 



/ M \, (Xi — Xj) + (yi — Vj) ■ The total potential energy of the system is 
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where the first sum in Eq. (j3J) is the potential energy of confinement, and the second sum 
is the potential energy due to particle-particle interactions. Equation (J3J) can be nondimen- 
sionalized to give 

where C/ is the characteristic potential energy, and & = Xj/r , r?j = j/i/r , and p^- = r.y/r 
are normalized distances. We define the characteristic length scale 

r o = — r:r- 6 

mo;^ 47re 

using the longitudinal oscillation frequency u 0x . The dimensionless parameters in Eq. (J5]) 
are the particle number n, the well anisotropy a 2 and the Debye shielding parameters k, 
where 



respectively. To compare this model to experiment, a 2 and k must be measured in the 



experiment. 

This model [Eq. (J5])] has three parameters: n, a 2 and k. The isotropic well is given by 
a 2 = 1. Without loss of generality, we assume that the anisotropic well has a 2 > 1 (i.e., 
uJq v > uJq x ) so that the major axis of the potential well lies in the x (longitudinal) direction. 
An unshielded Coulomb interaction corresponds to k — 0. As k increases the interparticle 
force becomes more localized. Given n, a 2 and k, a solution of the model [Eq. (jSJ)] is a set 
of particle positions that minimizes U. Properties of such solutions have previously 

been investigated by Candido, et al. [13]. To minimize U/U we use simulated annealing 



a. 



together with a final step of conjugate gradient minimization [7|. For a given configuration, 
normal modes and their associated frequencies can be computed from the dynamical matrix. 

A zigzag transition [ll] is a transition from a ID straight line configuration to a 2D 
configuration. For particle coordinates {xi,yi} measured with respect to the cluster's center 
of mass, the cluster's length and width can be characterized by the rms values 



n * — ' v n 

Consequently, a zigzag transition is a transition from y rms = to y rms > caused by a change 
in one of the model parameters. If a cluster is initially in a straight line configuration, then 
for constant a 2 and k a zigzag transition will occur as n is increased. We denote the critical 
value of n, which is the smallest value of n in the zigzag configuration, by n c . A 1D-2D 
transition also occurs when k is increased above a critical value k c . If a cluster is initially 

n 

in a 2D configuration, then increasing a 2 causes a transition to a ID cluster [13( for which 
Vi — Vi = above the critical value a 2 . As a consequence, ID configurations are independent 
of a 2 when a 2 > a 2 . 

An unbounded ID chain can be modeled by letting uoq x — > while ujQ y remains finite. 
Longitudina. confinement can be adneved either by nsing periodic bonndary conditions Q 
or a ring topology [22]. It is then convenient to define the characteristic length scale using 



the transverse frequency uo y 
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q2 (9) 



This gives a transverse Debye shielding parameter k t which is related to k [Eq. ([7j)] by 

Kt = t = m^ £k - (10) 

A zigzag transition occurs when the ID lattice constant a < a c where a c is a critical value 
19j, [22| . For an unbounded Yukawa chain, the dimensionless critical lattice constant a c /ror 



is a solution of [22] 



which depends only on the transverse shielding parameter kt- For a pure Coulomb inter- 
action, k t = 0, (a c /r 0T ) 3 = (7/4) ^(3) so that a c /r 0T ps 1.28. The critical lattice spacing 
decreases as the Debye shielding parameter increases (see Fig. [5] below). 
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Figure 1: Schematic of the experimental setup. Nearly identical spherical dust particles are confined 
in a biharmonic potential well created in the rectangular depression between four conducting bars 
placed on the rf powered electrode. Experiments were performed for confinement geometries with 
bar separations d = 25.4 mm and d = 14.0 mm. 

III. EXPERIMENT 



Dusty plasma experiments were performed in the Dusty Ohio Northern University experi- 



menT (DONUT) 



m 
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241 ] . An argon plasma was created around an 89-mm diameter 



powered electrode in a radio frequency discharge at 13.56 MHz. A blocking capacitor allows 
the electrode to develop a negative dc self-bias that levitates the negatively- charged dust 
particles. As shown in Fig. [TJ the biharmonic well is formed at the minimum of a confining 
geometry consisting of four rectangular aluminum bars placed on the powered electrode. The 
end bars measure 6.35 mm x 12.7 mm x 76.2 mm, while the inner bars measure 6.35 mm 
x 12.7 mm x 50.8 mm. The distance d between the two inner bars can be changed to vary 
the dimensions of the confining rectangular depression, and thereby change the anisotropy 
parameter a 2 . Clusters were made using monodisperse melamine formaldehyde spheres with 
a nominal diameter of 9.62 ± 0.09 /mi. As explained previously {9}, we believe that the dust 
particle diameter is closer to 8.94 ± 0.18 /im. 
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To determine dust particle positions, the particles are illuminated by a red diode laser 
and viewed using a 2/3 inch CMOS camera with a telecentric lens mounted above the top 
face of the electrode. For these experiments, we recorded 4097 frames of video at ~ 30 
frames/s for each particle configuration to determine center-of-mass (cm.) and breathing 
frequencies. A side-view camera was used to verify that out-of-plane motion was minimal. 

Two different confinement geometries were studied for similar plasma conditions. In the 
first, the inner bars were separated by d — 25.4 mm, while in the second they were separated 
by d = 14.0 mm. For the 25.4 mm x 50.8 mm well, the neutral Ar pressure was 12.4 mtorr 
(1.65 Pa), the rf power was ~ 10 W forward, the dc self bias on the electrode was —89.0 
V, and particle positions were recorded with a resolution of 16.51 /zm/pixel. For the 14.0 
mm x 50.8 mm well, the neutral Ar pressure was 12.1 mtorr (1.61 Pa), the rf power was 
~ 9 W forward, the dc self bias was —83.0 V, and positions were recorded with a resolution 
of 16.77 /im/pixel. 

Normal mode frequencies were determined by projecting the particle's thermal motion 
onto the center-of-mass and longitudinal breathing modes 10, l24j . A Fourier transform of 
the time history of the mode amplitude gives the power spectral density for that mode, 
which is that of a driven damped harmonic oscillator. For the neutral pressures used, 
the oscillations are underdamped and the power spectra display a clear resonance peak. 
Measuring the center-of-mass frequencies u 0x and u 0y directly determines the anisotropy 
parameter a 2 , while comparing the longitudinal breathing frequency for ID configurations 
to model solutions determines the Debye shielding parameter k. 



IV. EXPERIMENTAL RESULTS 



For the 25.4 mm x 50.8 mm confining well, nine sets of particles were analyzed with 
n = 2 to 19. Representative configurations are shown in Figs. [2(a)-(g). For n < 5 the 
particles are in a ID linear configuration. When one more particle is added (n = 6) the 
cluster changes to a 2D zigzag configuration, so n c = 6. As n increases the number of zigzags 
also increases until zigzags stretch from one end of the cluster to the other (n = 8, 9). For 
n — 19 the system displays a full (5,14) elliptical shell structure {13I 

For the narrower confining well, d = 14.0 mm, we analyzed twelve sets of particles for 
nine different values of n (2 < n < 17). Measured configurations are shown in Figs. Efa)-(g). 
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Figure 2: (a)-(g) Measured particle positions for confining bar separation d = 25.4 mm. (h)-(n) 
Computed positions for a 2 = 9.24 and k = 3. By matching y rms between the experiment and 
model, we find ro ~ 1.40 mm. Agreement between measured and computed configurations is 
excellent. For n < 5 the configurations are linear, at n = 6 a zigzag develops and for n = 19 a fully 
elliptical cluster with a well-defined shell structure is seen. Both measured and computed figures 
have a 1:1 aspect ratio. 

In comparison to d = 25.4 mm, we expect the anisotropy parameter to be larger so that 
the critical value of n is increased. For these conditions, clusters with n < 9 are in a ID 
configuration. A zigzag configuration is seen for n = 10, so n c = 10. As n further increases 
the zigzag region expands away from the center of the cluster. However, even for n — 17 
the cluster still has short linear tails at each end [13] and is not an elliptical configuration. 

The anisotropy parameter a 2 for each confining well was determined from measurements 
of the center-of-mass frequencies excited by thermal noise [h]] for the x and y directions, 
as shown in Fig. |H For the 25.4 mm x 50.8 mm well [Fig. H(a)], mode temperatures were 
found to be 300 — 400 K, indicating that the clusters are stable and in equilibrium with the 
neutral gas component. The cm. frequencies do not depend on the number of particles, 
so the clusters do not perturb the potential well. By averaging over the measured cm. 
frequencies, we find u 0x = 7.00 ± 0.06 rad/s and ujq v = 21.2 ± 0.1 rad/s. Equation (jBJ) then 
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Figure 3: (a)-(g) Measured particle positions for a confining well with d = 14.0 mm. (h)-(n) 
Computed positions for a 2 = 30.7 and k = 4. We estimate ro ~ 1.65 mm by comparing y rm s f° r 
the model and experiment. For n < 9 the configurations are linear, at n = 10 a zigzag develops. 
For n = 17 the cluster remains in a zigzag configuration. Both measured and computed figures 
have a 1:1 aspect ratio. 

gives a 2 = 9.24 ± 0.2, so that the anisotropy parameter has been precisely determined. 

For the 14.0 mm x 50.8 mm well [Fig. H](b)], data were taken at n = 6 for two different 
sets of particles and n = 9 with three different sets to estimate the spread in the measured 
frequencies, as can be seen in Fig. 11(b) . In comparison to d — 25 A mm, we observed a larger 
range of mode temperatures, 300 — 500 K, indicating that the clusters are somewhat less 
stable. The cm. frequencies are again found to be independent of n, and the average cm. 
frequencies are u> 0x = 5.37±0.06 rad/s and u 0y = 29.7±0.1 rad/s, giving a 2 = 30.7±0.7. In 
comparison to the d = 25.4 mm case, uo y has increased as expected, while u>o x has decreased 
slightly, even though the long side of the rectangular well (50.8 mm) has not changed. This 
indicates that decreasing d is pushing the sheath out of the concave depression formed by 
the bars. 

The Debye shielding parameter k was estimated by comparing measured longitudinal 
breathing frequencies for several of the linear configurations to normal mode frequencies 
calculated using the model. Since the breathing oscillation varies the interparticle spacing, 
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Figure 4: Measured center-of-mass frequencies and longitudinal breathing frequencies vs particle 
number n determined from thermally excited oscillations in (a) the 25.4 mm x 50.8 mm confining 
well and (b) the 14.0 mm x 50.8 mm well. Broken lines are average values, and the solid line in 
(b) is a linear fit to u>b r . 

it probes the dependence of the interparticle potential on particle separation, and therefore k. 
For an unshielded Coulomb interaction k — 0, the squared normalized breathing frequency 
(uJbr/uJox) 2 = 3 irrespective of n, and the unshielded Coulomb regime is k < 0.2. For the 
25.4 mm x 50.8 mm well with n — 4, the experimental value uj^ = 14.08 rad/s, so that 
(ufo I ujqx) 2 = 4.05, giving k « 2.6, while for n = 5, the experimental value u;& r = 14.27 rad/s, 
so that (ubr/^ox) 2 = 4.16, giving k ~ 3.1. The measured values of (uji, r /uj 0x ) 2 are clearly 
not consistent with k — 0, and we conclude that Debye shielding cannot be neglected when 
modeling these clusters. The uncertainty in k is fairly large, so we take k = 3.0. Using 
this value of k, we compare the measured cluster width y rms to the dimensionless cluster 
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width to find tq = 1.40 mm, q = —1.3 x 10 4 e and A = 0.47 mm. These values of q and A 
are consistent with measurements made in isotropic wells for similar discharge conditions 
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241 ]. In this case, the particle separation, which is a = 0.72 mm at the center of 



the n = 5 cluster, is greater than the Debye length, emphasizing the importance of Debye 
shielding. 

For d = 14.0 mm the Debye shielding parameter was estimated from the normalized 
breathing frequencies for clusters with n = 6, 8 and 9 particles. For n = 6 we find 
(u br /uj 0x ) 2 = 4.32 and 4.38, for n = 8, (00^ / u Qx ) 2 = 4.20 and for n — 9, (u br /uj 0x ) 2 = 3.98, 
4.24 and 4.78. From this data we estimate k ~ 4, giving r$ = 1.65 mm, q = —1.3 x 10 4 e and 
A = 0.41 mm. As we show in the next section, k = 4 is very close to the critical value for 
the zigzag transition, which may somewhat explain the spread in the breathing frequencies 
for n — 9. The physical parameters q and A are consistent with the values found for the 
d = 25.4 mm well even though k is somewhat larger due to the decrease in u 0x [Eq. (J7|)]. 

Equilibrium configurations computed from the model [Eq. (jSJ)] for a 2 = 9.24 and k = 3 
are shown in Figs. 121(h)- (n) for comparison to the experimental configurations. For each value 
of n, the experimental and predicted positions are very similar, and the particle arrangements 
are identical. In particular, the zigzag transition occurs at n = 6 in both cases, so that 
experimentally the critical value n c = 6. For a 2 = 9.24 and k = 0, n c = 7, which does not 
agree with the experimental results. A comparison between the measured configurations and 
configurations computed for a 2 = 30.7 and k = 4 is shown in Fig. [3l Again, the measured 
and computed configurations show excellent agreement and the particle arrangements are 
identical. Further, the critical value of n for the zigzag transition is the same for the 
experimental and the model results. For a 2 = 30.7 and k = 0, the critical value n c = 13. In 
the experiment we find n c = 10, so that our results are not consistent with k = 0. The very 
good agreement between experiment and model, and the consistency of the results for two 
different confining wells, indicates that the experimental results are robust. 

In Fig. [5] we compare the values of the experimentally measured lattice constant a for 
the last straight configurations, n = n c — 1, to the unbounded theory of Eq. ( TTTf) for the 
critical lattice constant. For both potential wells, the last straight configuration has n odd, 
so we approximate a by the average of the distances between the central particle and its two 
nearest neighbors. For d = 25.4 mm with n = 5, we find a = 0.72 mm, and for d = 14.0 mm 
with n = 9, we find a = 0.60 mm. Using the measured values of tq, k and a 2 we then 
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Figure 5: Critical lattice parameter a c /roT for an unbounded straight chain [Eq. vs the 

transverse Debye shielding parameter kt- The data points are the experimentally measured values 
for the last straight configuration for d = 25.4 mm (circle), and d = 14.0 mm (diamond). The 
measured points lie close to, but above, the stability curve in the stable region. 

calculate for d = 25.4 mm: a/r^T = 1-08 and kt = 1.43, and for d = 14.0 mm: a/r^T = 1-13 
and kt = 1.28. The experimental points lie close to the instability line, but slightly above 
it, in the stable region. Since for both cases a/r T < 1.28, decreasing kt (e.g., increasing 
A) while holding a/r T constant moves the cluster into the unstable region, causing a zigzag 
transition. Here finite size effects do not appear to be very important, which may be because 
for kt > 1 the zigzag instability is dominated by nearest neighbor interactions 22]. 

In Fig. [6] we compare the measured cluster width y rm s with model solutions as a function 
of n for the measured values of a 2 and k. The only adjustable parameter is the length 
scale for the cluster r , which was chosen to give good agreement between the model and 
experiment. In both cases we see an abrupt increase in the cluster width which is associated 
with the zigzag transition. Above the transition, the data exhibit a power law behavior 
which is consistent with a continuous phase transition. Agreement between the model and 
the experiment is quite good. For d = 25.4 mm there is a second structural transition at 
n = 15 which corresponds to the change from a zigzag configuration to an elliptical shell 



configuration [Fig. [2(g)] 
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This transition may be roughly analogous to the transition 
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Figure 6: Dependence of cluster width y rm s on particle number n comparing experiment and model, 
(a) Experimental data for bar separation d = 25.4 mm and model solutions with a 2 = 9.24 and 
k = 3 scaled using tq = 1.65 mm. (b) Experimental data for bar separation d = 14.0 mm and 
model solutions with a 2 = 30.7 and k = 4 scaled using tq = 1.40 mm. The dashed lines are power 
law fits to the model points. 

from two to three parallel chains in the unbounded system Such a transition is not 

seen for d = 14.0 mm since n is not large enough, as confirmed by the fact that the n = 17 
configuration is an extended zigzag [Fig. EJ^g)]. 



V. DIMENSIONAL PHASE TRANSITIONS 

A phase transition is a sudden change in some property of a system, called an order 
parameter, due to a small change in a control parameter. Within this conceptual framework, 
the zigzag transition in these clusters can be viewed as a dimensional phase transition 
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We characterize the cluster size in 



between one- dimensional and two-dimensional states 
the longitudinal (x) and transverse (y) directions by the rms values of the particle positions in 
the respective directions [Eq. flSJ)]. In particular, y rms is a good choice for an (unnormalized) 
order parameter since y rm s — in the ID configuration and y rms > in the 2D (zigzag) 
configuration. The state variables that determine the system configuration are then n, k 
and a 2 , where n is discrete and k and a 2 are continuous. 

Figure [7| demonstrates that n c decreases as k increases using a 2 = 30.7. That is, expres- 
sions which predict the critical value of a 2 for a given n with k = [13y are incorrect when 
the strength of Debye shielding is such that the interaction is not essentially unshielded. 
Interestingly, even though n is small and discrete, the cluster width above the transition is 



well characterized by a power law 



19] 



cx(n-<r (12) 



where n' c is a continuous critical n, and n c = \n' c ] . That n' c is continuous indicates there may 
be a continuum theory for the zigzag transition where n is also continuous. For a 2 = 30.7 
and k = 0,1 and 4, we fit the first five points after the transition to find n' c = 12.34, 10.82 and 
8.99, and a critical exponent v = 0.430, 0.387 and 0.310, respectively. Here v decreases with 
increasing k. As discussed above, k = 4 is very close to the critical value when a 2 = 30.7, 
and we find n' c very close to an integer value. In fact, for k = 4, n c = 10 so we expect n' c > 9, 
which is not satisfied here due to a small uncertainty in the fitting coefficients. 

The computed dependence of cluster length and width for n = 5 and 6 and a 2 = 9.2 on 
the Debye shielding parameter k is shown in Fig. [HI For the finite model, in contrast to the 
unbounded case [19j, increasing k decreases the nearest neighbor distance, and therefore the 
linear particle density. For n = 5 the critical value of k for the zigzag transition is k c = 4.22. 
That is, the cluster is in a ID configuration for k < 4.22. For k > 4.22, y rms is positive 
and increases rapidly with k. The cluster length x rms has a discontinuous first derivative at 
the phase transition. The inset of Fig. [HJ shows the transition for n = 6, where the critical 
value is k c = 0.45. Consequently, for n = 6 and k < 0.45 the cluster is linear. Since we 
find experimentally that the n = 6 cluster is in the zigzag configuration, we conclude that 
the experimental Debye shielding parameter must lie in the interval 0.45 < k < 4.22. For 
k > k c , the cluster width has a power law behavior 

Vrms OC (k /? c ) , (13) 
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Figure 7: Cluster width y rm s/ r Q vs particle number n for Debye shielding parameters k = 0, 1 
and 4 with well anisotropy a 2 = 30.7. The critical value of n at which the zigzag transition occurs 
decreases as k increases. Solid lines are power law curves fitted to the first five points following 
the zigzag transition. 

where (3 is a critical exponent that is independent of the normalization of y rrns - For the 
n = 5 and 6 cases illustrated in Fig. [SJ we find (3 = 0.463 and 0.450, respectively. This 
analysis was repeated for a 2 = 30.7 with n = 9 and 10. For n = 9, k c = 4.08 and (3 = 0.463, 
while for n = 10, k c = 1.77 and j3 = 0.469. Experimentally, we find n c = 10, which means 
that for the experiment k must lie in the interval 1.77 < k < 4.08. The critical exponent for 
the zigzag transition vs k is nearly the same for the four cases considered here, so it may 
be that in Yukawa clusters there is a universal critical exponent (3 ~ 0.46 for the zigzag 
transition caused by changing the Debye shielding parameter k. 

For fixed values of n and k, a 2D-1D transition (an "inverse zigzag") takes place as the 
anisotropy parameter a 2 increases, as shown in Fig. [9] for n = 5 and k = 3. As a 2 increases, 
the dimensionless cluster length increases while the width decreases. Near the transition, 
the width exhibits a power-law approach to y rrns = 0. Assuming 



we find the critical value a 2 = 8.74 and the critical exponent 7 = 0.387. That is, if k — 3 
and a 2 < 8.74 then n = 5 particles will be in a zigzag configuration. These results are 
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Debye shielding parameter k 

Figure 8: Computed cluster width and length for anisotropy parameter a 2 = 9.2 for n = 5 and 
(inset) n = 6 particles vs the Debye shielding parameter k. Solid lines show the power law fit to 
Eq. (|13|) . while the dashed line has been added to guide the eye. For n = 5 particles there is 
a critical value k c = 4.22 below which the cluster is one dimensional and above which it is two 
dimensional. For n = 6, k c = 0.45. 

consistent with the experiment where we observed a ID configuration for a 2 = 9.24. When 
a 2 > a 2 the model results are independent of a 2 , which can be seen by the constancy of 
x rms , since yi — r)i — and the configuration only depends on n and n. 

VI. CONCLUSIONS 

We have studied one- and two-dimensional Yukawa clusters with a small number of par- 
ticles n < 19 confined in biharmonic potential wells both experimentally and theoretically. 
Experiments were performed in the Dusty Ohio Northern University experimenT (DONUT). 
For n less than a critical value n c , the clusters are in a one-dimensional straight line state. 
When n = n c the cluster undergoes a zigzag transformation to a two-dimensional state. In 
our experiments, the anisotropy of the confining potential well was accurately determined 
by measuring the frequencies of center-of-mass oscillations excited by thermal noise in both 
the x (longitudinal) and y (transverse) directions and the Debye shielding parameter was 
estimated from the measured longitudinal breathing frequency of ID clusters. Experimental 
and model data show excellent quantitative agreement, confirming that dusty plasma is a 
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Figure 9: Computed cluster width and length for n = 5 with shielding parameter k = 3 vs 
anisotropy parameter a 2 . A 2D-1D phase transition is with a critical value a 2 = 8.74. The fitted 
power law (solid line) gives a critical exponent 7 = 0.387. 

very good real-world system for studying ID and 2D Yukawa systems and the transitions 
between these states. In particular, strongly-coupled linear configurations with n < n c are 
true ID systems for which normal modes are either purely longitudinal or purely transverse. 

Our results clearly show that Debye shielding is important for our experimental condi- 
tions, and our results are not consistent with physics in the unshielded Coulomb regime. For 
a given value of the potential well anisotropy, the critical particle number n c decreases as 
the shielding parameter k increases (i.e., as the Debye length decreases). For the measured 
well anisotropies, the experimental values of n c are below those predicted for an unshielded 
Coulomb interaction (k = 0), indicating k > 0. This is reinforced by noting that the mea- 
sured Debye lengths are less than the particle separation and that the normalized squared 
longitudinal breathing frequencies (uj br / uj 0x ) 2 > 3. 



2l| for a similar experiment 



Our finding that k > contradicts the conclusion in Ref. 
where the results where said to be consistent with an unshielded Coulomb interaction. The 
method used in Ref. 



2JJ to determine cluster parameters is a static analysis that treats both 
the anisotropy parameter a 2 and shielding parameter k as free parameters and compares the 
observed n c with that predicted for k = 0. When the decrease in n c with k described in the 
present work is considered, it seems likely that static analysis method 21] is only weakly 
constrained and cannot be used to accurately determine cluster parameters. 
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Finally, we have demonstrated that the zigzag transition in a Yukawa cluster can be 
viewed as a phase transition from a one-dimensional state to a two-dimensional state. 
Though this was previously demonstrated for unbounded systems jl9|. here the number 
of particles is finite and really quite small. This is true for transitions initiated by changing 
the Debye shielding parameter, the potential well anisotropy and the number of particles. In 
all three cases, we find that the transverse cluster width has a power law dependence near the 
transition, indicated that transition behaves as a continuous phase transition with a critical 
exponent. For 1D-2D transitions caused by increasing the Debye shielding parameter, we 
tentatively identify a universal critical exponent [Eq. ( fl3j) ] with a value (3 ~ 0.46. 
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